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Abstract. I give a brief review and interpretation of the evolution, mass 
loss and pulsation of AGB stars in the Large Magellanic Cloud. 



1. The Magehanic Clouds in an astro-physical context 

Our view of the place of mankind in the Universe changed dramatically with the 
recognition by Abbe (1967) that the Large and Small Magellanic Clouds (LMC 
& SMC) are in fact nearby galaxies composed of stars and nebular matter, 
separate from our own Milky Way galaxy. Four decades later, the discovery 
in 1907 by Henrietta Leavitt of a relation between the pulsation period and 
apparent brightness of Cepheid-type variables in the LMC counts amongst the 
most important discoveries in modern astronomy. 

Over the past few decades, a great deal of astro-physics has been learnt 
from studies of stars in the Magellanic Clouds, for several reasons: (i) stars 
in each of the Clouds are all at the same distance to us, and therefore their 
relative luminosities are well determined; (ii) this distance is well known; (iii) 
the Clouds are near enough to study individual stars; (iv) interstellar extinction 
to and within the Clouds is small; (v) as dwarf galaxies, the Clouds are fairly 
big, and they have experienced star formation over much of their history. Hence, 
plenty of stars can be found in different evolutionary stages, including stars in 
clusters for which accurate ages may be determined, (vi) The LMC and SMC 
have metallicities of ^ and ^ solar, allowing the study of metallicity effects. 
However, as the analyses and astrophysical problems they address become more 
detailed, these advantages meet their limitations: (i) the depth of the Clouds 
has been measured and is not always negligible; (ii) different methods yield 
a distance modulus of the LMC ranging from ~ 18.1 to 18.7 mag (1 adopt 
distance moduli to the LMC and SMC of 18.5 and 19.0 mag, respectively); (iii) 
current instrumentation does not always provide the required spatial resolution 
or detection sensitivity; (iv) although extinction is much less than in the plane of 
the Milky Way, it is not always negligible even at infrared (IR) wavelengths, and 
can be very patchy in places; (v) the number of stars in very rare, short-lived 
evolutionary phases is limited — sometimes to less than one at any moment 
— hampering statistical studies. Also, the star formation history of the LMC 
is complex, with a lack of 4 to 10 Gyr-old stellar clusters. And (vi) there are 
no (identified) samples of either metal-rich stars, with [Fe/H]> 0, or extremely 
metal-poor stars, with [Fe/H]<C —2. 
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2. Asymptotic Giant Branch (AGB) stars 

The AGB phase is the final stadium in the life of a star with a Zero- Age Main- 
Sequence (ZAMS) mass of ~ 0.5 to 8 Mq. The lower limit to this mass range 
is set by the requirement that the temperature and pressure at the bottom of 
the hydrogen mantle surrounding the electron-degenerate core (oxygen and car- 
bon/neon) is sufficiently high to ignite hydrogen shell burning, driving the evo- 
lution of the star along the asymptotic giant branch towards higher luminosity 
(growing core) and lower surface temperature (expanding mantle). The AGB 
evolution terminates either because the core has reached the Chandrasekhar 
mass limit of Mmax ~ 1.4 Mq beyond which the core becomes gravitationally 
unstable, the corresponding luminosity of -Lmax ~ 55, 000 Lq defining the tip of 
the AGB, or because the mantle has diminished to the extent that the hydrogen- 
burning shell becomes exposed and quickly runs out of fuel, leaving the hot and 
compact core (white dwarf) to cool into oblivion. The upper limit to the mass 
range of AGB progenitors is set by the condition for the ignition of carbon 
burning in the non-degenerate core of supergiants. 

The hydrogen-burning shell produces an underlying helium shell which ig- 
nites as soon as sufficient pressure has built up, thereby causing the mantle to 
expand and the hydrogen burning to switch off. Such Thermal Pulses (TP) 
re-occur on timescales from At ~ 10^ yr for the most massive AGB stars up to 
At ~ 10^ yr for the least massive AGB stars. The TP itself is brief (~ 1% of At), 
but it takes long (up to 30% of At) for a star to recover from the luminosity 
dip of up to a factor two whilst the hydrogen shell burning gradually takes over 
from the helium shell burning. The effect of the TP on the mantle and emergent 
luminosity is relatively stronger for less massive stars. The TP- AGB phase lasts 
for t =fewxl0^ yr, up to t =fewxlO^ yr for the most massive AGB stars unless 
mass loss is severe. 

Energy transport from the nuclear burning site to the stellar surface of an 
AGB star proceeds via convective motion through the mantle: a rising gas ele- 
ment continues to rise. Much theoretical understanding needs to be gained about 
the exact criterium for convective instability and the description of convective 
fiows and boundaries, especially because convection not only transports energy 
but also matter and hence affects the stratification of elemental abundances 
throughout the mantle. When during a TP the convection zone penetrates into 
layers containing the products of nuclear burning, these heavy elements will 
become mixed into and become over-abundant in the stellar photosphere where 
light leaves the star. This is the third time in the evolution of intermediate-mass 
stars that such dredge-up event may happen, and probably occurs only for the 
more massive AGB stars (M ^ 1.5 Mq). Products of the slow neutron-capture 
process, s-process elements such as zirconium, technetium and barium are used 
to diagnose the TP-AGB history of stars. Carbon may become enhanced too, 
greatly affecting the chemistry in the cool stellar atmosphere (Teg ^ 4000 K) 
if after binding carbon and oxygen into carbon-monoxide (CO) not oxygen but 
carbon remains to form more complex molecules. Such carbon stars arc easily 
distinguished from oxygen-rich AGB stars as their spectra are dominated by 
absorption bands of CN, HCN, C2 and C2H2 rather than TiO, VO and SiO. 
However, nuclear burning may take place at the inner boundary of the convec- 
tion zone, this Hot Bottom Burning (HBB) converting carbon into oxygen and 
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nitrogen and preventing the formation of a carbon star. HBB occurs in massive 
AGB stars (M ^ 4 Mq) and enhances their luminosity, but switches off once 
the mantle has nearly gone due to severe mass loss. 

The stability criterium for motion within the mantle not only determines the 
occurrence of convection (turbulent motion) , but also the occurrence of pulsation 
(organised motion). In the hydrogen ionization zone not far below the stellar 
surface, and in the first helium ionization zone, the absorption coefficient of the 
gas to radiation from the stellar interior is very large; hence the internal energy 
becomes high and the mantle expands until the opacity drops and the mantle 
shrinks again. These radial pulsations levitate the atmosphere high enough for 
dust formation to occur and radiation pressure to blow away the circumstellar 
dust — taking with it the collisionally coupled gas. Mass-loss rates of M ~ 10~^ 
to 10~^ M0 yr~^ cause the pre-mature termination of AGB evolution, exposing 
the white dwarf after the star has shed its mantle. 

3. AGB stars in the Magellanic Clouds 

Searches for AGB stars in the Magellanic Clouds have followed three main strate- 
gies: (i) monitoring surveys to find Long- Period Variables (LPVs), mostly of 
Mira type defined by regular pulsation on the order of a year with a visual am- 
plitude A\r > 2.5 mag; (ii) spectroscopic or photometric surveys to find bright 
M-type or carbon stars, and (iii) mid-IR surveys to find AGB stars with circum- 
stellar dust emission. 

3.1. Optically bright AGB stars 

Amongst the most influential systematic searches for magellanic M-type giants 
and carbon stars are the study of selected fields in the bar of the LMC and in 
the SMC (Blanco et al. 1980); Blanco & McCarthy 1983; Progel & Richer 1983), 
the construction of the carbon star luminosity function in the LMC (Costa &; 
Progel 1996), and the survey of magellanic clusters for which ages and hence the 
ZAMS masses of the cluster members can be determined (Aaronson &; Mould 
1985; Progel et al. 1990; Westerlund et al. 1991). 

Although much of the physics remains in paramcterised form — notably 
convection and dredge-up — evolutionary models can be tuned to reproduce 
the carbon star luminosity function quite well (Marigo et al. 1999). When these 
results are confronted with the luminosities and chemistry of cluster AGB stars in 
the Magellanic Clouds (Progel ct al. 1990), the match is not entirely satisfactory 
(Pig. 1 left) due to both imperfections in the models as well as difficulties in 
deducing the observed quantities. 

The first grand survey for variable stars resulted in the catalogue of Harvard 
Variables (Payne-Gaposchkin 1971) which, however, were mostly supergiants. 
Subsequent surveys identified pulsating AGB stars in sufficient numbers for a 
P-L relation to be proposed (Glass &; Lloyd-Evans 1981) and refined (Wood 
et al. 1983; Peast et al. 1989; Hughes k Wood 1990). Mira's did indeed seem 
to obey a remarkably tight and unique P-L relation (Fig. 1 right) suggesting a 
robust instability threshold and a single pulsation mode, and their possible use as 
distance indicators, but already a few complications may be noticed: (i) The P-L 
relation has a different slope for periods shorter/longer than P' ^ 400 d, which 
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Figure 1. Left: models (Marigo et al. 1999) and cluster data (Frogel 
et al. 1990); Right: P-L relation in the LMC (Feast et al. 1989; Hughes 
& Wood 1990) and evolutionary tracks (Vassiliadis & Wood 1993). 



evolutionary tracks (Vassiliadis & Wood 1993) link to a ZAMS mass M' ~ 2 
Mq. Previous evolution, for instance whether (M < M') or not (M > M') the 
star experienced violent ignition of core-helium burning at the tip of the first 
red giant branch (RGB), may be the cause for a difference in internal structure 
and pulsational properties. Also, (ii) carbon stars and especially red stars have 
somewhat larger periods for their luminosities, which may be due to their larger 
radii as they have suffered significant mass loss, (iii) A paucity of luminous AGB 
stars, in particular carbon stars with L ^ 20, 000 Lq corresponding to M 4 
Mq, was explained as due to HBB (Iben 1981; Wood et al. 1983). This was 
confirmed by the identification of lithium-enriched S-type stars (C/0 just below 
unity) with 20, 000 < L < 55, 000 L© (Smith et al. 1995). 

3.2. Dust-enshrouded AGB stars 

Studies of optically bright AGB stars are biased against stars in the final phases 
of AGB evolution that have become dust-enshrouded as a result of mass loss. 
It is difficult to study their properties as their stellar surface is obscured at 
optical wavelengths and the photosphere at IR wavelengths becomes dominated 
by thermal emission from the circumstellar dust. 

The IRAS mission in the early 1980s provided a catalogue with a few hun- 
dred mid-IR point sources in the Magellanic Clouds that are candidate dust 
shells around AGB stars (Loup et al. 1997). Follow-up on such IRAS sources 
and similar dust-enshrouded supergiants (Wood et al. 1992; van Loon et al. 
1998, and references therein) culminated in the detailed quantitative study of 
~ 50 dust-enshrouded AGB stars in the LMC observed with ISO (van Loon et 
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Figure 2. P-L diagram in the LMC (left), and compared with the 
SMC and the galactic bulge, centre and globular clusters (right). 



al. 1999). The main results are that (i) luminous carbon stars do exist amongst 
dust-enshrouded AGB stars, confirming the HBB scenario with mass loss; (ii) 
the threshold above which HBB occurs is confirmed to be 4 M© (van Loon 
et al. 2001); (iii) winds of AGB stars are dust-driven (e.g. Gail & Sedlmayr 
1986); (iv) mass- loss rates exceed nuclear burning mass consumption rates and 
reach a maximum Mj^^xOcL in a distinct superwind phase; (iv) mass-loss rates 
increase with lower Tefj — this effect is strongest during the early stages of cir- 
cumstellar dust formation — but may be only little dependent on metallicity for 
[Fe/H] ^ — 1 (van Loon 2000). Great improvement in our knowledge of AGB 
mass loss is expected from recent surveys of the Magellanic Clouds at near-IR 
(2MASS, DENIS) and mid-IR (MSX, ISO) wavelengths, and their follow-up. 

3.3. Large-scale variability surveys 

Much progress is made in detecting variables in the LMC thanks to micro-lensing 
surveys (MACHO, EROS, OGLE). The P-L diagram (Fig. 2, left) not only shows 
the classical Mir a sequence but at least three more roughly parallel sequences. 
Two of these, at shorter periods and populated by semi-regulars (many of which 
are regular but have smaller amplitudes than Miras), are believed to correspond 
to overtone modes with the Miras pulsating in the fundamental mode (Wood & 
Sebo 1996; Wood et al. 1999). This interpretation is supported by the behaviour 
of pulsation amplitudes &: velocities (Scholz & Wood 2000), period & radius, and 
non-linear models (Ya'ari &; Tuchman 1996, 1999) for galactic Miras, but if the 
mixing length describing convection is allowed to vary then Miras could still 
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pulsate in the first overtone (Barthes 1998). Another sequence at longer periods 
than the Mira P-L is explained as binarity (Wood et al. 1999). 

It is estimated that not all but ^ 65% of the TP-AGB stars are LPV, of 
which (only) ~ 10% are Miras (Cioni et al. 2001), and the evolution in the P-L 
plane is complex: the pulsation is a function of mass and time — not necessarily 
in a continuous or even monotonous fashion. It is then remarkable how consistent 
the P-L diagram is populated by stars in the LMC, SMC and galactic centre, 
and the galactic bulge (M <^ 2 Mq) and globular clusters (M ^ 1 Mq) despite 
an order of magnitude difference in metallicity (Fig. 2, right). 

An overdensity of stars at the two short-period sequences just below the 
maximum luminosity at the RGB indicates these may be RGB stars, whilst 
above the maximum luminosity at the AGB, supergiants are found on sequences 
similar to but slightly displaced from the Mira P-L and the long-period sequence. 
This suggests that certain similarities in the internal structure of RGB stars, 
AGB stars and (red) supergiants cause similarities in their pulsation. 

Carbon stars populate the Mira P-L where stars with Mzams ~ 1-5 to 4 
M0 are found, as expected from third dredge-up and HBB. Carbon stars on the 
sequence immediately at shorter periods than the Mira P-L, at L ~ 8,000 Lq, 
may have Mzams ~ 3 to 4 Mq, but the continuation of the evolutionary track 
is devoid of stars until it reaches the Mira P-L at a twice as high luminosity. 
Perhaps, after having become carbon stars on the Mira P-L, these stars are 
temporarily semi-regulars, possibly due to a TP. Surprisingly, the most luminous 
carbon stars obey the Mira P-L. Thought to have become carbon-enriched after 
an extended period of severe mass loss terminated HBB, they would be expected 
to have the longest periods of all — longer than on the Mira P-L. 



4. The connection betw^een pulsation and mass loss of AGB stars 

When AGB stars arrive on the Mira P-L, their mass-loss rates have become so 
severe that further evolution is along tracks of constant luminosity, with the 
radius (and hence pulsation period) growing as the structure of the mantle re- 
adjusts itself to the decrease in mass. Hence the dust-enshrouded stars (Wood 
1998; Whitelock et al., in preparation) are on the Mira P-L or to longer peri- 
ods. Describing the Mira as a polytrope {R\fM = constant) and the pulsation 
as a harmonic oscillator (27r/P ~ 1/ifrccfall — 9{R)/R)-, the period evolves as 
Pocl/M. For instance, a Mira with P = 250 d (Mzams = 1.5 Mq) would end 
up on the long-period sequence at P ~ 750 d and M ^ 0.5 Mq, which is a 
reasonable white dwarf mass. The longest observed period for an evolved equiv- 
alent of a Mira with P = 700 d (Mzams = 5 Mq), however, has P ~ 1400 d and 
M ~ 2.5 Mq. Such star must still lose more than a solar mass, but stars with 
correspondingly long periods (P ~ 3000 d) have not (yet) been found. 

The growth in radius/period of mass-losing AGB stars may become trun- 
cated when the dynamical timescale, P, approaches the thermal timescale, 
tc^GM'^/LR, and the mantle no longer pulsates adiabatically. This happens 

when Pc^^irM/aT^g, corresponding to P ~ 2000 d for a star of M = 2.5 Mq 
and Tcfj = 2500 K (note the moderate dependence on stellar parameters). By 
this time, the surface gravity is so low that mass loss could be initiated without 
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Figure 3. K-band amplitudes and pulsation energy input rate, versus 
period (left) and luminosity (right) in the LMC, SMC and Milky Way. 



the help of pulsation: the star may quite literally be boiling over, and this could 
be responsible for the superwind phase required to shed the remainder of the 
mantle that may become visible as a planetary nebula. The ability of the mantle 
to restore thermo-dynamic equilibrium on a dynamical timescale may also be 
relevant for the faint part of the long-period sequence (L ^ 10,000 Lq). Some 
disturbance, for instance tidal interaction in a binary, could trigger instabilities 
on a thermal (not orbital) timescale, causing irregular mass loss that results in 
obscuration events, super-imposed on a "traditional" (adiabatic) pulsation in 
the fundamental mode (Mira P-L) or overtone mode (shorter-period sequences). 

K-band amplitudes (a measure for the bolometric amplitude) generally in- 
crease with period, but for P ^ 400 d the increase becomes steeper (Fig. 3, left). 
This could be related to the different slope of the Mira P-L relation for periods 
shorter /longer than P' ~ 400 d. The diagram suggests that, for longer periods, 
carbon stars and stars in the Magellanic Clouds pulsate less vigorously than 
galactic oxygen-rich stars. However, the pulsation of a luminous star with small 
amplitude in magnitude may well involve more energy than that of a faint star 
with large amplitude. The time-averaged rate at which radiative energy is stored 
in the mechanical pulsation, E = L/2x (exp[Ai^/2.5] — l)/(exp[AEr/2.5]-|-l) has 
a theoretical maximum of E^i^i/ L = 0.5: the maximum energy rate available 
during that half of the pulsation cycle in which the mantle is lifted is L. The 
observed E increases with luminosity (Fig. 3, right), saturating at E'max/-^ = 0.4 
— i.e. very efficient — with no clear distinction between stars of different chem- 
istry. As the maximum mass-loss rate was found to increase proportionally to 
luminosity (van Loon et al. 1999), I obtain Mmax oc Emax- This explains the 
absence of a clear metallicity-dependence of the mass- loss rate (van Loon 2000). 
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